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I. INTRODUCTION 


One problem that affects computations in quantum 
field theory (QFT), in both its applications to condensed 
matter or high energy physics, is the possible appear¬ 
ance of infrared (IR) divergences in the massless or crit¬ 
ical cases [T]. In this class of theories, the breakdown of 
the mean field, or Landau, approximation was one of the 
original problems that triggered the development of the 
renormalization group (RG) theory j^. This lead to the 
concept of theory space, running couplings and ultimately 
the understanding of universality. We then learned how 
to compute universal quantities like critical exponents or 
universal ratios but, ultimately, a technique free of IR 
divergencies that allows the computation of correlation 
functions, which are the main object of interest in QFT, 
did not emerge. Only later conformal field theory (CFT) 
was developed to fill this gap, but only at criticality, i.e. 
at fixed points [3. Also, applications of CFT techniques 
have been limited to two dimensions and only recently 
progress has been made in three dimensions, although it 
has been mostly numerical [ 1 ]. 

In this paper we want to fill this gap by showing how 
to compute correlation functions, or more precisely one- 
particle-irreducible (IPI) vertices (generated by the ef¬ 
fective action), along a RG trajectory, in such a way that 
if this trajectory is critical or massless and hits a fixed 
point in the IR, no divergencies occur and a perfectly fi¬ 
nite analytical result emerges. If instead, the trajectory 
leads to a massive theory in the IR, then still the compu¬ 
tation goes through. This reflects the fact that a theory 
is indeed a trajectory in theory space that connects a mi¬ 
croscopic, or bare action, to the effective action in the IR, 
which contains the physics. 

The idea is to perform the path integral by weighting 
the momentum modes that contribute to it according to 
their RG relevance, i.e. we weight each mode according 
to the value of the running couplings at that scale. In 
other words, we apply the RG improvement procedure 
also to the correlation functions by integrating their RG 
flow with the prior knowledge of the phase space struc¬ 


ture, i.e. with the running couplings as known functions 
of the scale, functions obtained by integrating their re¬ 
spective beta functions. After explaining the method in 
abstract, we will show how it concretely works in the case 
of a model with non-trivial phase diagram as the three 
dimensional Z2-scalar theory. 


II. THE METHOD 


The fundamental reason for the success of RG com¬ 
putations resides in the RG improvement, which usually, 
and originally, is applied only to the couplings: here we 
will apply it to the full effective action. The starting point 
it the trace-log formula for the one-loop effective action: 

T = S+^TilogS‘^^'>, ( 1 ) 


where S is the bare action and 5'G) jg Jtg Hessian. The 
trace on the r.h.s. stands (in fiat spacetime) for an inte¬ 
gration over coordinate space, or momentum space, and 
for an integration over any representation of internal and 
spacetime symmetries that the fields may carry. Obvi¬ 
ously eq. has to be regularized to make sense: let’s 
use a simple mass regularisation /r for the IR and a cutoff 
regularisation A for the ultraviolet (UV). In this way the 
effective action will become dependent on /r through the 
mass regulator, while the bare action and the trace will 
implicitly depend on A: 


= -S'a -I- ^TrA log 




( 2 ) 


Note that all the A dependence arising from the trace 
is compensated by counter-terms contained in the bare 
action S'a and so the effective action depends only on /i. 

Now let’s pause a moment to remember how we usu¬ 
ally implement the RG in the case of a coupling, say A 
dimensionless. We reabsorb the UV divergencies into the 
renormalized coupling: = Aa — GA^^ log ^. To obtain 

the beta function we then take a derivative with respect 
to to obtain = CX\, then we perform the RG 





improvement on the r.h.s. Aa —>■ A^ to finally obtain 
the RG equations M^^A^ = CAJJ. This is basically the 
original procedure of Gell-Man and Low 

It is at this point that we will deviate from the standard 
practice and introduce our method by repeating the RG 
procedure but at the level of eq. (i). For convenience, 
from here onwards we change notation and we will denote 
the RG scale with k instead of /r. By taking a scale 
derivative of eq. (§ we get 

1 

where t = logfc/fco and fco is some reference scale. At 
this point we perform the RG iinprovement substituting 
—>■ on the r.h.s. of eq. ^ 


III. THE MODEL 

We start by considering real scalar theories in d- 
dimensional Euclidean space with Z 2 ~symmetry defined 
by the following bare action: 

Sa = J . (7) 

We are interested in computing the effective action Fq 
using eq. Q and eq. as described in the previous 
section. The effective action is a complicated non-local 
action; here we will perform an expansion in the field ip 
and retain all terms up to order ip'^. We thus fix the form 
of the scale dependent effective action Tfc to be 


dtTk 



dtRk 
+ Rk ’ 


(4) 


where we stepped ahead and introduced the cutoff kernel 
function Rk = In general, this kind of cutoff is not 
sufficient to ensure the trace to be finite. In order to make 
eq. not UV divergent, Rk is usually taken to be a func¬ 
tion of the Laplacian and required to fast fall-off for large 
values of its argument. In this case also its t-derivative 
decays fast for large argument values and this makes the 
trace in the r.h.s. of eq. 0 convergent. In eq. @ the 
form of Sk is the same of the bare action Sa with all bare 
couplings replaced by their running counterparts. 

Knowing the flow of Ffc, in analogy to the case of the 
couplings, it is possible to compute the renormalized ef¬ 
fective action by integrating eq. Q between the UV and 
IR 


ro = rA 


1 f^dk dtRk 
2 io k^^s^^^+Rk’ 


(5) 


where Ta = ^a is the bare action which is taken as the 
initial condition of the flow in the UV and Tq = T/j is 
the renormalized effective action. In order not to create 
confusion, we want to stress that in this formalism the 
subscripts "0" naturally denotes renormalized quantities 
and will be often traded for the subscript "i?". 

The action T/j is the generating functional of IPI ver¬ 
tices 


{ip{xi)ip{x2) ■ ..ip{Xn)) IPI 


6ip{xi) . ..6ip{Xn) 


. ( 6 ) 

ip—0 


Finally, one can do an even further step and switch on 
all possible renormalized couplings, that we imagine to 
be equal to zero in the bare action, and make the sub¬ 
stitution —>■ in eq. in this way obtaining a 

closed equation for the running effective action F^. One 
can show that this equation is exact [S]. 


Tk = jd’^x + Sfc(-n) -I ml]ip+^ip'^Fk{-n)ip‘^'^ , 

( 8 ) 

where Efc(—□) is the running polarization function, 
is the running mass of the scalar field and Ffe(—□) is the 
running four-point structure function. Here □ = 9^9^ 
stands for the fiat space Laplacian. The constant part 
of the four-point structure function defines the running 
quartic scalar coupling A^, i.e. Ffc(O) = Afc. The full 
effective action in principle contains also all terms with 
a higher field power but, since they are not relevant for 
our discussion, we neglect them and restrict just to the 
truncation of eq. (§. 

We now take the Hessian of the bare action in eq. Q 
and make the RG improvement of the couplings: tua —t 
nik and Aa —t Afc. We find: 


=-D + ml + ^p>^, (9) 

which is the quantity that has to be inserted in eq. Q. 
In order to properly account for threshold effects, it is 
convenient to choose the argument of the cutoff Rk{z) to 
be z = —□ -\-^ip^. In this way eq. 0 assumes the form: 


dtTk = \T^hkUu + ^p>\ml 


( 10 ) 


where we defined the threshold function hk{z,uj) = 
dtRk(z) ^ need now to evaluate the functional trace 


z+uj+R k(z) 

in eq. (lOl; this can be easily done using the non-local 


heat kernel expansion, which we review in the Appendix. 
The use of the expansion reported in eq. (A4| leads us 
to: 










dtVk = - 




ds hkis, mf.) Tr e 


1 1 


2 (47r)^/2 
1 1 
2 (47r)^/2 


J d‘^x J dshk{s,m1)s 


1 _ ^s(^2 hAs^ip‘^fij(sX)(p‘^ + 0{ip^) 


Afc 


Qd[hk]-—Qd_j^[hk]^+^ip diQd_ 


^xai-0 


+ 0{ip^) 


( 11 ) 


where X = — □. In the last line we have used the ex¬ 
plicit expression of fu given in eq. ( |A5[ ) and we intro¬ 
duced the Q-functional notation, also described in the 
Appendix. For computational convenience we have also 
dropped, on the r.h.s. of equation (11), terms arising 


from the scale derivative dtRk present in the function hk 
and proportional to the beta function of \k- The neglec- 
tion of these terms renders the flow equation (111 similar 


to the proper-time RG equation studied, for example, in 
|7]. For a detailed discussion about the implications of 
this approximation see |H|- 

Plugging in the l.h.s. of eq. 0 the ansatz for Tk 
given in eq. ([^ and comparing the same held monomials 
on both sides, we can read off the flow equations: 


dt^kix) = 0 


dtm^. = -- 


1 Afc 


2 


dtFkiX) = I- 


2 (47r)^/2 


Q d_^ \hk\ 


d^ Q d_2 




,( 12 ) 


where we defined h^(z,a;) = hk{z + a,uj). The last line 
of eq. (12) reduces to the flow equation of the quartic 
coupling Afc in the limit A —> 0: 


dtXk = - 


^ Qd-2 [^k] ■ 


2 (47r)‘^/2 


(13) 


The four-point structure function is computed integrat¬ 
ing the flow of Fk{X) in eq. (12) from k = A (UV) down 
to fc = 0 (IR): 


Wt- 

AF{X) = ^dtFk{X) 

(14) 

The function Fo{X) = Fa{X) — AF{X), where Fa(A) = 
Aa is the UV initial condition, completely determines the 
IPI four-point correlation function: 

{(p{xi)ip{X2)ip{x3)(p{x4)) IPI 

= l:d^xi5xx2Fo{-D)6xx3Sxxi + permutations; 

D 

connected correlations can then be reconstructed via 
standard relations. 



Eq. (141 is our main result so far, and before moving 


on to consider an explicit application, we will stop a mo¬ 
ment to explain its meaning and in particular how it is 
related to the procedure we announced in the introduc¬ 
tion. Consider the massless case mfc = 0 for simplicity 
and let’s undo the RG improvement of the self-interaction 
coupling for a moment by restoring Aa- Using the fact 
that hf. = dt \og{z+a+Rk) = dtl% we can rewrite eq. (14) 
as: 


where 


AF{X) = J ^Aa i9t I one-loop I (15) 


one-loop = 


2 (47r)^/2 


d^ d _ 2 




(16) 


is the standard one-loop result in presence of the cutoff 
Rk, as the reader can easily convince itself. Now we can 


explain the crucial point of our construction: in eq. (16) 
one can perform the RG improvement Aa —>■ Afc after 
exchanging the A:-integral with the bare coupling, or be¬ 
fore and then being unable to do the switch. In the first 
case we simply reproduce the standard one-loop compu¬ 
tation since formally = 1, while in the second 

case we perform an RG improved computation, here of 
the structure function Fq{X), where we effectively weight 
each momentum shell dk with the value of the running 
coupling at that point, here A^. This second procedure 
is the one we propose and its application is encoded in 
eq. (141. As we will now show with an explicit exam¬ 
ple, this procedure is more powerful, consistent and well 
defined than the standard perturbative approach since it 
automatically encodes the information about of the phase 
space trajectory we are integrating along, through the 
scale dependent couplings, here Afc and rrik for which we 
have previously solved. 


IV. THE PHASE DIAGRAM 


We restrict ourselves to scalar theories in three dimen¬ 
sion, which are UV finite but suffer of IR divergencies in 
the massless-critical case and are perfect to show how our 
method works. The beta functions for and Afc can be 




























The values of the two fixed points are slightly shifted: 



Figure 1: Phase diagram in the three dimensional case. One 
can see the Gaussian and Wilson-Fisher fixed points and the 
trajectory connecting them, along which we integrate the flow 
of the structure function F’fe(a;). 


read off by setting d = 3 in eq. (|12| and eq. (|13|: 


1 


dtm^ = -X 


2 (47r) 


Q1 [^fc] — 






Q \2 o\2 

dtXk = 


(17) 


2 (47r)5 Stt^ + mf 

We have used the explicit values of the Q-functionals 


given in eq. (B 61 and eq. (B71 of the Appendix. The beta 
function for the dimensionless couplings mj. = m\/k‘^ and 


Afe = Afe/fc are: 


II 

0™2 Afe 1 

^ d-TT^ 1 - 1 - ml 


3X1 1 

dtXk = 


There exist two fixed 

points at: 

= -A 

A* = (IR) 

ml = 0 

A* = 0 (UV) 


(18) 


The IR fixed point is the well-known Wilson-Fisher fixed 
point and has one attractive and one repulsive direction. 
The UV fixed point is the Gaussian one and has both 
directions attractive, as shown in Fig. 1. The critical ex¬ 
ponents at the IR fixed point can be easily computed, see 
for example [1]. It is possible to find an analytic solution 
for the flow if one considers eq. (181 in the approxima¬ 


tion in which the mass in the denominator on the r.h.s. 
is neglected: 


dtml = 
dtXk = 


- 2 wfc - 


Afe 

47r^ 


t 3A2 


(19) 


~ 2 1 
< = -3 

* 

II 

loo 

1 

(IR) 

ml = 0 

0 

II 

* 

(UV) 


but their properties are not. The analytic solution of the 
system in eq. (191 is: 


Afe = 


C 2 2Ci 

k 

SClTT^ 

3Ci+k ’ 


&cl 


log(3C'i -I- k) 


( 20 ) 


where Ci and C 2 are integration constants with, respec¬ 
tively, dimension of mass and mass squared. Notice that 
the UV fixed point is reached, for k —>■ 00 , independently 
of the values of these constant, since it is attractive in 
both directions. On the other hand, the IR fixed point is 
reached, for fc —>■ 0 , only if a particular relation between 
the integration constants holds, i.e. only if we are on the 
critical trajectory. 

We want to repeat here the fundamental point: a the¬ 
ory is represented by a trajectory in theory space. All 
trajectories can be parametrized by the value of the renor¬ 
malized mass 


mji = lim k^mj. = C 2 + log(3C'i). 

k—>-0 


( 21 ) 


The unique trajectory that flows between the two fixed 
points is defined by: 


m% = 0. 


( 22 ) 


Using this condition in eq. (20 1 gives the following solu¬ 
tions: 


o 2Ci f3Ci + k 

= -^ + ^log 




Afe — 


SClTT^ 

3Ci + k 


(23) 


Eq. (231 describes the unique trajectory in the {fh\,Xk) 


plane that flows from the Gaussian fixed point in the 
UV to the Wilson-Fisher fixed point in the IR. In Fig. 1 
this finite trajectory is shown in red. All non critical 
trajectories are characterized by the condition 7 ^ 0 . 
Positive values of label trajectories that flow towards 
the symmetric phase, while negative values of label 
trajectories that flow towards the broken phase. 

We will now solve the theory along the critical line 
connecting the Gaussian and Wilson-Fisher fixed points 
and compare the result with the standard perturbative 
treatment. 


V. THE FOUR-POINT FUNCTION 


The four-point structure function is obtained by eval¬ 
uating the expression in eq. (141 for d = 3. We have 


AF{X) 


3 1 f°°dk^2 k^ fX\ 

^TxJo T ^k^+mf[¥) 


(24) 























where we have used the explicit form of the Q-functional 


given in eq. (B8l to compute 


g{u) = log -y/u) + log - 2). 

(25) 

This is the analytic result we obtained by using the op¬ 
timized cut-off in eq. (B41. Other cutoffs may be imple¬ 


mented and in general will give different forms for g{u), 
i.e. g{u) is a regulator dependent quantity. We also re¬ 
moved the UV cutoff A —>■ 00 since the theory is finite in 
this limit. 



Figure 2: The dashed curve is the standard perturbative re¬ 
sult {X, Xr,0) in the massless limit, the dotted curve 

represents the improved structure function AF^{X, Xr/Stv^) 
in the mk = 0 case and the solid curve represents the full im¬ 
proved structure function AF^\X, Xr/Stt^ ,0). All the three 
structure functions are evaluated for Xr = 4/%/3. 


A. One—loop computation 


The structure function at one-loop is obtained by 
performing the RG improvement only after having ex¬ 
changed the bare coupling and the scale integral in 
eq. (24 \ and by setting the mass to its renormalized value: 






1 


87 r 2 Jo 


dk 

/ 

^xx 

lo k 

fc2 -I- 



3X1 , ^[X 

-^Arctan--, 


(26) 


where vriR and Xr are the renormalized couplings. In the 
massless limit eq. (261 reduces to: 




16 Va' 


(27) 


Eq. (26) and eq. (27|are known to be the correct per¬ 
turbative results mm. Notice that in the perturbative 
case the result is not sensitive the the RG trajectory we 
are on, i.e. there is no way to codify this information 
and the result is the same in all cases. In particular, the 
massless (critical) expression in eq. (27) is IR divergent 
at large scales X —)■ 0, signalling the breakdown of 
Landau-mean field theory. In general, for d < 4, the IR 
divergences at itir = 0 are more and more important as 
the perturbative order increases. If p = ^/X represents 
the momentum scale, the four-point function, computed 
at one-loop through a "bubble" diagram, behaves as 
p-{4-d) n times inside UV convergent 

diagrams, will produce IR divergencies if n(4 — d) > d. 
For more details see [T]. We will see in the next section 
that the RG improved computation of the four-point 
function structure function (bubble diagram) will cure 
this problem. 


B. RG improved computation 


We now perform the RG improved calculation where 
the scale integration contains the running couplings 
and Afc evaluated on the critical trajectory connecting 
the two fixed points. Every momentum shell will thus be 
weighted accordingly to the strength of the interaction 
and, differently from the perturbative case, we will obtain 
an IR finite four-point structure function. 


For a better analytical understanding, we will first in¬ 
tegrate the flow of the structure function eq. (24) in the 
case nik = 0, which we will call AF^. In this case theory 
space is one dimensional and the complete fc-dependence 
of Xk is given in eq. (23). We then find the following 
result: 


AF\X,Ci) = 


r k 

VX Jo (3Ci + fc)2 




2 , ASCiVX , 6Ci _ Vx + 6Ci 




2 , 


Vx + 6Ci 
VX-6Ci 


. 4PolyLog 


6Ci -f VX 
6 Ci - VX 


(28) 


The value of the integration constant Ci can be fixed by imposing that AF^ {X, Ci) has to reproduce the pertur- 




































bative result in the UV limit X ^ oo: 


3.0 


AF^{X,Ci) 


12Tr‘^Cf 


JC —>-<50 


(29) 


Comparing with eq. (271 we then have: 


12 - ^ 


c. = ^. (30) 


In the IR limit X —0 the improved structure function is 
given by: 



X 


{X,X]i/8Tr^) = Xr 


(31) 


and is finite as should be. The improved structure func¬ 
tion AF^ {X, Afl(/87r^) differs from the standard one-loop 
perturbative result AF^^{X,X[i,0) of eq. (271 at small 
values of X. A comparison between the perturbative and 
the improved computation is shown in Fig. 2, where we 
have set Xr = The improved result takes correctly 

into account the fact that the theory flows to an IR fixed 
point and thus it furnishes a IR finite quantity. In the 
language of perturbation theory this is equivalent to say 
that the bubble diagram is finite and no more divergent 
as XjxfX in the IR. This is the exemplification of how the 
RG cures the breakdown of Landau-mean field theory if 
properly implemented, as does the method we propose. 

We now relax the condition rufc = 0 and we integrate 
the flow of the structure function in eq. (241 substituting 
for m\ and their complete fc-dependence as given by 
eq. (pOl. We have: 


Ai^^^(A,Ci,m^) = 


87r2 Jo 


dk k 


XI , 

1 + ^fc^l 


f X 


where C 2 has been traded for using the relation in 
eq. (|21| ). We first perform a numerical integration of 
eq. ( [ 32 I along the critical trajectory (m|. = 0) setting 
Cl = Xr/8'k‘^ = l/-\/l27r^ in order to compare the re¬ 
sult directly with the previous computations, as shown in 
Fig. 2. Again the four-point structure function is finite 
in the IR, namely for A —>■ 0. The integration of eq. (321 
along non critical trajectories with > 0 provide also a 
finite IR limit of the structure function. Fig. 3 shows the 
behavior of the structure function obtained by numeri¬ 
cal integration of eq. (321 considering different positive 
values of m\ and setting Ci = I/vT^tt^. Trajectories 
that flow to the broken phase (to|j < 0 ) make the inte¬ 
grand of eq. (321 to develop a singularity along the inte¬ 
gration path and we can not describe this situation prop¬ 
erly within our truncation. In this case the truncation 
we use turns out to be inconsistent because it describes 
a theory which is quantized around the wrong vacuum. 
At the fixed point the theory is a CFT, so the A —0 
limit can be compared with what is known from CFT. 
Consistently we find that the IPI part of the four-point 
function, i.e. with the four propagator stripped-off, is to 
first approximation constant El- 


Figure 3: The plot represents different behaviors of the 
structure function AF^^(A, l/\/T27r^, m|{) computed, start¬ 
ing from below, for mfj = 1, 0.1, 0.01, 0.001 and 0. 


VI. OUTLOOK 

In this paper we showed that it is possible to extend 
the standard RG improvement of the couplings to an im¬ 
provement of the full effective action by generalising the 
original procedure Gell-Mann and Low. This method can 
be used to make any standard one-loop computation in¬ 
clude the non-trivial information encoded in the phase di¬ 
agram of the theory. In this way it is possible to cure both 
UV and IR divergencies when appropriate fixed points are 
available as, respectively, starting or ending point of the 
RG flow. We also showed how this makes one-loop com¬ 
putations become sensitive to the RG trajectory along 
which one integrates the theory, as it should. 

We exemplified the method with an application to 
a three dimensional real scalar with Z 2 -symmetry and 
showed how to obtain a consistent IR finite four-point 
structure function as the result of the integration along 
the trajectory connecting the Gaussian (UV) fixed point 
with the Wilson-Fisher (IR) fixed point. 

Our method can be applied to any QFT, for example to 
improve a mean field analysis. It also works the other way 
around, i.e. going to the UV; in this case one can obtain 
finite results in theories with an UV fixed point. For ex¬ 
ample, in four dimensions one can apply these ideas to the 
computation of the gravitational effective action along a 
trajectory controlled by an UV fixed point; a first step in 
this direction has been done in [12]. Other applications 
can be to NLcrM m or to controllable asymptotically 
safe theories El- 
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Appendix A: non—local heat kernel 


we can rewrite eq. (Al \ as: 


We now show how to compute the functional trace in 
eq. (flOl: 


dt^k = ^Tihki-O + U,oj), 


(Al) 


with U = Introducing the Laplace transform 

hk{s,uj), 


hk{z,uj) = / dshk{s,oj)e 

Jo 


dtVk = 


ds hk{s, a;)Tr e 


-s(-n+c/) 


(A3) 


The trace in the r.h.s of eq. (A31 is computed using the 


non-local heat kernel expansion m which retains the 
infinite number of heat kernel coefficients in terms of non- 
(A2) local structure functions or form-factors and reads: 

I 


Trg-s(-n+c/) ^ J d^x^l^l-sU + UM-sa)U + 0{U^)Y 


(A4) 


where the structure function fu{x) is given by: 

= (A5) 

Appendix B: Q—functionals 

The Q-functionals are defined as 

pOQ 

Qn[K]= dss-^e-^‘^hk{s,uj), (Bl) 

Jo 

where h‘^{z^uj) = hk{z + a,uj). For a positive real number 
n > 0 it is possible to rewrite the Q-functional as a Mellin 
transform m- 


Qn{hk] — 


1 


r(n) 


dzz" '^hk{z,uj), (B2) 


while for a negative half-integer n = — < 0 we have 

m- 

C? 2 tti + 1 




We compute the Q-functional using an optimized cutoff 
of the form 


R^[z) = (fc 2 - z)e{k^ - z ). 


(B4) 


In this case the function hk{z, uj) takes the following form, 

2Pe{k^ - z) 


hk(z,uj) = 


z + LU + (k^ — z)0(k^ — z) 


(B5) 


We have 


1 1 

Q^i^k] = / dzz~^hk{z,uj) 

1 ( 9 ) Jo 


r(i)jo 

4 




(B 6 ) 


poo Jm+i 

(B3) and 


^ poo' id 

Q_i [hfc] =- -j= / dz z~^ —hk{z + a,u}) 

^ V'^ Jo dz 

2 9[kJ — a) 

VJr + w - a 


(B7) 


d^Q_i 




2 P f\, e{k^-xai-0) 
Vlf fc 2 + ml Jo v'fc 2 - A^(l - ^') 

2 P 1 


■P 


Vx 


log -X) + log ^ + ^S (A - 


2k-VX 


Vx - 2 k 


(B 8 ) 


This completes the evaluation of the Q-functionals. 
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